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ABSTRACT 

Radial velocity searches for extrasolar planets have discovered stars undergoing slow 
accelerations. The accelerations are likely due to one or more gas giant planets or 
brown dwarfs orbiting with period longer than the total time span of observations. 
The stellar acceleration is proportional to the mass of the companion divided by 
the square of its radius from the star. In this paper we predict the distribution of 
accelerations using a Monte Carlo method and assuming a semi-major axis and mass 
distribution for the companions. Radial velocity surveys suggest that 20 and 10% of 
stars surveyed exhibit accelerations above ~ 10 and 25 m s~^ yr~^, respectively. We 
find that an extrapolation of the size and period distribution found by radial velocity 
surveys within ~ 5 AU or companion imaging surveys predicts too few systems with 
detectable accelerations. The fractions of stars with accelerations above these values 
can be matched if every star hosts a gas giant planet or brown dwarf with semi-major 
axes between 4 and 20 AU or if a significant fraction (~ 70%) of stars host 3 or more 
giant planets in this range. If brown dwarfs companions are not discovered in imaging 
surveys then we would infer that there is a large population of gas giant planets present 
with semi-major axes in the 5-30 AU range. 



1 INTRODUCTION 



Radial velocity searches for extrasolar planets often 
find linear trends in the observed veloc it y of a star 



jButler et all l2006l: Iwittenmver et al 



l2007l : iPatel et all l2007l : iMarcy et al.l boOSi ). These 



20061: Iwright et all 



corre- 
spond to velocity variations or constant accelerations of the 
star by massive objects in the system. As the time span 
of ob servations increases these acceleration s begin to vary 
(e.g.. IWright et al.ll2007l : iPatel et ai]|2007D . then finally a 
substantial fraction of the an orbital period becomes rec- 
ognizable as the entire Keplerian velocity curve is recog- 
nized. Trends or accelerations h ave been seen in radial ve- 
locity surveys (|Butler et al ]|2006l ) and used to place limits on 
mEiss and radii of long period planets residing in th e system 
jWittenmver et al.ll2006l. l2007l : IWright et allboOTl ) . Planets 
are not considered detected by radial velocity or transit sur- 
veys until at least a full orbital period is seen l|Q'Toole et al.l 
|2008| ). This is justifiable as it is difficult, perhaps impossible 
when noise is present, to constrain simultaneously planet 
mass, semi-major axis and orbital eccentricity from radial 
veloci ty observations alone over only a fraction of a pe- 
rio d (iBrownl I2OO4I: IWright et all bOOTl ). but see the study 
bv lPatel et al.l (|2007| ). Once an entire period is observed, 
the period of the orbit sets the planet's semi- major axis. As 
the amplitude of the radial velocity perturbation depends 
on both the mass and period, once the period is known, 
the mass of the planet can be estimated. Because only poor 
constraints on both planet and semi-major axis can be made 



from constant or slowly changing accelerations, radial veloc- 
ity surveys place only poor constraints on exoplanets at large 
semi-major axis where the orbital periods are long. 



In the next decade we have the prospect of both more 
sensitive radial velocity surveys (nearing Im/s in preci- 
sion) and multi-object radial velocity surveys. These surveys 
will increasingly find stars with measured accelerations that 
don't correspond to complete orbital periods. In this paper 
we consider the distribution of these trends and how they 
can be used to constrain the mass and semi-major axis dis- 
tribution of massive objects in outer exo-planetary systems 
even though the exact masses and semi-major axes are not 
directly measured. A similar approach can b e adopted fo r 
astrometric surveys for long-period planets (|011ing| [20071 ). 
By outer we mean above 5 AU. The orbital period at 5 AU 
around a solar mass star is about 11 years, corresponding 
to the lifespan of a decade length multi-object radial veloc- 
ity survey. This semi-major axis is within the limit of most 
current direct planet detection methods. However the statis- 
tics of objects inferred from the statistics of radial velocity 
trends will be complimentary to and motivate direct detec- 
tion surveys sear ching for brown dwarfs and planets at wider 
separ ations (e.g.. lMetchev fc Hillenbrandll2008l : iKraus et al.l 
|2008| ) and transit surveys that primarily find planets within 
an AU of their host star. 
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2 THE ACCELERATION OF A STAR CAUSED 
BY A LOW MASS COMPANION 

We first consider tiie acceleration of a star caused by a single 
planet or brown dwarf in orbit. In the center of mass frame 
of 2 gravitationally bound bodies with masses mi and m2 at 
positions xi and X2, the position of the more massive body 
is 



Xl 



m2 



X2 



(1) 



where M = mi + m2 is the sum of the masses of the two 
bodies, the mass ratio fi — m2/M and the vector between 
the two bodies f — xi — X2. The acceleration of the more 
massive body or host star 



Xl = fir. 

For the two body Kepler problem 
GMf 



(2) 



(3) 



where G is the gravitational constant and r is the unit vec- 
tor between the two masses. Combining the previous two 
equations we find that the acceleration of the central star 

= (4) 

The acceleration toward the viewer measured from changes 
in the radial velocity is ar — xi ■ eios where eios is the unit 
vector pointing toward the observer. In physical units 



187m s 



" y"^" (xu) 



(iF^) 



(5) 



We note that the magnitude of the acceleration depends on 
fj,M/r^ — Mp/r'^, where Mp is the mass of the planet and is 
independent of the stellar mass. We expect the distribution 
of accelerations to depend on the ratio of the companion 
mass divided by the square its semi-major axis. This depen- 
dence is reflected in the u pper limit plots (mass vs period) by 
IWittenmver et al] l|2006l ) (see their Figure 4) derived from 
trends (accelerations) seen with the McDonald Observatory 
planet search program. For comparison the amplitude of the 
radial velocit y var iation during an entire orbit is propor - 
tional to /i-^/M/a (e.g., equation 1 ICumming et al.1l2008h . 
We expect the statistics of accelerations to differ from those 
of radial velocity amplitudes. 

The distribution of extrasolar planets discovered in ra- 
dial velocity surve ys is consistent with dN/dloga oc a*^ * 
l|Marcv et ahllioosl ) implying that there is increasing num- 
bers of giant planets at larger semi-major axes. For a semi- 
major axis greater than 10 AU, where an entire orbital pe- 
riod would not be detected in a moderate time span ra- 
dial velocity survey, equation [S] predicts an acceleration of 
1.9m s~^ yT~^{Mp/Mj) which is near the detection limit 
of radial velocity surveys. Here Mj is the mass of Jupiter. 
However companion m asses in the brown d warf range (be- 
tween 13 and 80 Mj, iBurrows et al.l [l997h at semi-major 
axes between 10 and 100 AU can cause accelerations de- 
tectable in radial velocity surveys. Hence accelerations mea- 
sured from radial velocity surveys are sensitive to massive 
planets and brown dwarfs in the outer regions of extrasolar 
systems. This range of masses contains brown dwarfs and so 
the region identifled as the br own dwarf desert where ther e 
is a deficit of known objects l|Grether fc Lineweaveill2006h . 



Previous studies have recognized that incomplete orbits de- 
tected in radial velocity s tudies are most likely caused by 
compani on brown dwarfs ( Patel et al.ll2007l ) . 

, Gumming et all l|200^ ') found that ~ 23% of the to- 
tal stars surveyed showed a significant slope or velocity 
variations (corresponding to accelerations) across 8 years 
of observations. The mean error of their survey was less 
than 3 m/s so the velocity variations detected would be 
similar to accelerations (along line of sight) are abo ve few 
m s^^ yr^^- With less precision iNidever et all l|2002l ) found 
that 12% of the stars surveyed varied by more than above 
100 m s~^ across 4 years, corresponding to accelerations 
above 25 m s~^ yr~^- A period of 4 years corresponds to 
a semi-major axis of 2.5 AU. If we restrict the semi-major 
axis to greater than 3 AU and use Equatio n [51 the level of 
accelerations seen bv iNidever et al.1 (|2002l ') requires planet 
masses above 1 Mj if they are found at semi-major axes 
above 3 AU. We have conservatively chosen a semi-major 
axis corresponding to a period only slightly longer than the 
survey length. As the acceleration switches sign every half 
period it likely the objects responsible for these accelera- 
tions are at semi-major axis with periods larger than twice 
the survey length. For the longer Keck planet search 8 years 
corresponds to a semi-major axis of 4 AU. If we restrict the 
companion semi-major axis to 6 AU then again a planet 
mass above 1 Mj is required to account for the observed 
velocity variations. The large fractions of objects exhibiting 
detected accelerations suggests that there is significant pop- 
ulation of massive objects in the outskirts of exo-planetary 
systems. This is consistent wit h the extrapolat i on to larger 
semi-major axes discussed by I Gumming et al.l (j2008) who 
estimated the number of systems with planets based on size 
distribution of inner planets exhibiting full periods during 
the span of observation. 



3 MONTE CARLO SIMULATIONS 

To compare the number of objects with accelerations to 
those predicted from the size and period distributi ons we 
use a Monte Carlo method (similar to the study bv iBrownl 
l2004h . We compute the line of sight component of the accel- 
eration Or, based on a distribution of systems with planets. 
We first consider the distribution of accelerations that would 
be seen from a very simple population. We consider a group 
of systems, each with one massive planet or brown dwarf 
companion at random positions in their orbits and at ran- 
dom inclinations. Mean anomalies and arguments of perihe- 
lia are randomly chosen from a flat probability distribution. 
Eccentricities are chosen from a flat eccentricity distribution 
with minimum and maximum eccentricities emin and emax- 
Orbital inclinations are chosen by generating cos i with a 
flat probability distribution between -1 and 1 l|Brownll2004l . 
equation 17). After randomly generating orbital elements 
the Cartesian coordinates of each system are computed us- 
ing Kepler's equation. Then the coordinates of the system 
are projected into viewer coordinates using the randomly 
selected system inclination. Equation [5] is used to predict 
the line of sight component of the acceleration. Acceleration 
distributions are computed using 10^ sample star systems. 

Figure 1 shows the distribution of the radial or line of 
sight component of acceleration Ur computed for a popula- 
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tion of systems with 10 Mj bodies at 5 and 10 AU chosen 
with flat eccentricity distributions. For each semi-major axis, 
we show eccentricity distributions that consist of objects in 
circular orbits, range from to 0.5 or to 0.99. As planets on 
high eccentricity orbits have higher accelerations when they 
approach the star, the acceleration distributions extend to 
higher velocities when the eccentricity distribution is wider. 

For these simple systems we ask what is the fraction of 
objects that would be observed with accelerations greater 
than a particular value, 5, 10 or 25 m s~^ yr"^, /s, /lo and 
/25. These fr actions mimic the fracti o n of objects with trend s 
reported by ICumming et al.l (120081 ): iNidever et all (|2002l ). 
These fractions are listed in Table [1] and are measured from 
the same simulations that made the acceleration distribu- 
tions shown in Figure [T] 

The fraction of stars, /s, /lo and /25, with accelera- 
tions above 5, 10 or 25 m s^^ yr~^ is not high. The frac- 
tio n /25 does no t appr oach the 12% based on the survey 
bv lNidever et al.] l|2002l ) except for the systems with masses 
above 5 Mj mass planets at small radii at or below 5 AU. 
The fractions /s or /lo do not approach the 23% estimated 
(based on the report bv lCumming et al]|2008h unless masses 
are above 5 Mj. The fractions shown in Table [T] were com- 
puted assuming that every system contains a single massive 
planet. The fraction of stars showing accelerations is neces- 
sarily lower than the fraction of systems with planets (here 
assumed to be all systems) as some of the systems are viewed 
face on and systems are viewed at different times in their or- 
bit. The larger fraction of stars with accelerations compared 
to those predicted with this simple model suggests that one 
massive object is required at 5-20AU in almost all systems 
or more than one massive object is required in a significant 
fraction of systems. 

3.1 Accelerations predicted for a distribution of 
companion masses and periods 

We now consider Monte Carlo simulations for systems with 
planets and brown dwarfs of different masses and semi-major 
axes. The semi-major axis and planet mass distribution is as- 
sumed to be power law in form and is drawn from a specified 
range with limits amin, dmax and Mp^min and Mp^max- This 
type of distribution function was used to fit the radial veloc- 
ity planet search measurements at the Lick and K eck obser- 
vatori es and the Anglo Austra lian Telescope bv iMarcv et al.l 
(|2008l '): ICumming et al.1 (120081 ). We assume that the proba- 
bility of a planet with semi-major axis a and planet mass 



Mp 

dP = CM^a^ dM, 



da de. 



(6) 

l|Marcv et al.ll2008l ) find that a fa --1.1 and (5 ^ -0.6. The 
exponen t = —0 .6 is consistent with dN/dloga (x a'^'*. 
[Marcv et al.l (|2008l ) find a broad eccentricity distribution ex- 
terior to a = 0.1 AU with maximum eccentricity e = 0.93. 
The coefficient C is such that ~ 10.5% have planets with 
mass in rang e 0.3 to 10 MJ and orb i tal period in range of 2 
to 2000 days (|Cumming et al ] |2008l ). ICumming et al.l (|2008l ) 
computed the distribution for Mp sin i instead of Mp, how- 
ever o nly a small correction is needed as (sini) = n/A = 
0.785 (|Cumming et aL I I2OO8I . section 3.0). 

We now look at the distribution of accelerations for a 
population of a fixed companion mass but with a distri- 
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Figure 1. Probability distributions for the radial component of 
acceleration for a population of systems each with a 10 Mj planets 
at a semi-major axis of 5 AU (large points) and 10 AU (smaller 
points). The orbital inclinations, mean motions and arguments 
of perihelia are randomly chosen. Objects are drawn from 3 ec- 
centricity distributions for each planet semi-major axis. The first 
distribution is for circular orbits, the second with a flat distribu- 
tion ranging from eccentricity e = to 0.5 and the third randing 
from e = to 0.99. Higher eccentricities allow closer approaches 
and so higher accelerations. A lower mass planet would caused 
similar curves but shifted to the left because the accelerations 
would be lower. 



bution of semi-major axes. We consider randomly oriented 
systems populated by 10 Mj planets with semi-major axis 
chosen from a probability distribution dP oc a*^ , a flat eccen- 
tricity distribution ranging from e = to 0.99 and minimum 
and maximum semi-major axes between amin and amax ■ Fig- 
ure [5] shows that the probability distribution of the radial 
acceleration component is only weakly dependent on the ex- 
ponent of the semi-major axis distribution and on the inner 
radius of the distribution. These curves would move to the 
left for a population of lower mass planets as they would 
cause lower accelerations. 

We now consider a planet mass and semi-major axis 
distribution consistent with equation [S] We begin by nor- 
malizing the distribution to be consistent with the frac- 
tion of systems, ~ 10%, that have planets with mass be- 
tween 0.3 and 10 Mj and periods between 2 and 2000 days 
(|Cumming et aL I l2008l ). Each systems is assumed to have at 
most 1 companion. Figure [3] shows the acceleration distri- 
butions using this normalization. These probabilities now 
represent the likelihood of detecting systems with different 
radial components of acceleration. 

In Figure |3] we show 5 probability distributions. The 
comparison model is the c entral set of points on Figure [3] 
with a = -1.1, P = -0.6 (|Marcv et al.|]2008l ). This model 
has minimum and maximum companion semi-major axes of 
3 AU , amax — 100 AU. As shown from the curve just 
below the comparison model, the distribution is not sensitive 
to the minimum semi-major axes. This is likely because the 
distribution is normalized so that it extends at small semi- 
major axes to a probability consistent with the number of 
systems with planets discovered with radial velocity surveys. 
The comparison model has minimum and maximum planet 



4 



T3 




0.5 1 1.5 2 0.5 1 1.5 2 

log^Q acceleration (m s"^ yr"^) log.io acceleration (m s"^ yr"^) 



Figure 2. Probability distributions for the radial component of 
acceleration for a population of planets of a given mass but for 
different semi-major axis distributions. The comparison model, 
(solid circles) has 10 Mj mass planets chosen with a semi-major 
axis distribution dN oc a'' with /3 = —0.6. The probability for the 
comparison model is computed assuming that all systems have 
10 Mj planets with this distribution and with semi-major axes 
in the range between a^in = 3 AU and amax = 100 AU and 
eccentricities in the range to 0.99. The probability distribution 
is not sensitive to amax ■ Also shown are distributions for a higher 
minimum semi-major axis Qmin = 5 AU and shallower exponent 
fS = —0.4 with all other parameters the same as the comparison 
model. The probability distribution is only weakly sensitive to 
the inner radius of the distribution and the exponent of the un- 
derlying semi-major axis distribution. Power law distributions in 
the line of sight component of acceleration are predicted. 

masses of 0.1 and 100 Mj. The lowest curve shows the effect 
of truncating the mass distribution to planet masses below 
10 Mj. The upper two curves show the effect of changing 
the exponents a and /3. Shallower mass distributions have a 
larger number of more massive objects. This leads to a larger 
number of objects predicted with larger accelerations. 

3.2 Mass and semi-major axis likelihoods for a 
given acceleration 

For a measured acceleration such as that measured for GL 
849 with ar = 4.6 ± 0.8 m s'^ yi'^ (|Butler et all 120061 ') we 
can use our Monte Carlo code to measure the number of sys- 
tems sampled from a distribution that would be measured 
with this acceleration. The result of this using the distribu- 
tion with limits and exponents given in the first line of Table 
[2] is shown with contours in a plot of Mp vs a in Figure U 
Figures |3K a-nd b show contours for an eccentricity distri- 
bution with Cmax = 0.99, and 0.3, respectively. The highest 
probability occurs along a line proportional to Mp/a^ or 
Mpjc? « 0.03Mj/(AU)^ as e xpect ed from the scaling in 
equation [5] IWittenmver et aL iUustrated upper lim- 

its on planet mass as a function of semi-major axis with 
this scaling. Here we find that the most likely planet mass 
and semi-major axis mass also have the same scaling. The 
probability distribution is extremely broad implying that 
the mass and semi-major axis of the planet cannot be tightly 
constrained from the probabihty distribution alone. The dis- 
tribution is less broad if the eccentricity distribution is trun- 



Figure 3. Normalized probability distributions for the radial 
component of acceleration for a population of planets that varies 
in both mass and semi-major axis. The distributions have been 
adjusted so that they are consistent with the fraction (~ 10%) 
of systems with planets with mass between 0.3 and 10 M j and 
periods between 2 and 2000 days llCumming et aI.ll2008^ . Each 
system is assumed to have at most one companion. The central 
curve shows the comparison model that has parameters listed 
in the first line in Table [2] The other curves show probabilities 
when a single parameter is changed from those of the comparison 
model. Lowering of the exponent associated with semi-major axis 
(/3) raises the probability of detecting objects with low accelera- 
tions. This is because the distribution is constrained at high ve- 
locity accelerations by the radial velocity distributions. Lowering 
the exponent of the mass distribution (a) raises the probability 
of detecting objects at all accelerations as there are more mas- 
sive planets. Truncating the planet mass distribution lowers the 
probability of detecting all accelerations in the range shown here. 



cated. If the eccentricity distribution is broad then 10-100 
Mj brown dwarfs between 15-30 AU are somewhat more 
probable than gas giant planets. However, HST/NICMOS 
observations of this object in 1997 would have detected a 
brown dwarfs at a large separation, suggesting that the ac- 
celeration is caused by a gas giant planet closer to the star. 

3.3 Fractions of objects detected with specific 
acceleration values 

In Table [2] we list the fraction of systems with radial compo- 
nent, of their accelerations above 5, 10 and 25 m s~^ yr"^. 
Table [2] shows that these fractions are not strongly depen- 
dent on the inner semi-major axis of the distribution, amin so 
they are not strongly dependent on the time span of the ra- 
dial velocity survey. The fractions are weakly dependent on 
the lower planet mass Mp^min and the outer semi-major axis 
a-max of the distribution. However these fractions are depen- 
dent on the outer mass of the distribution and halve if the 
distribution is truncated at 10 Mj instead at 100 Mj. The 
exponent, a, setting the distribution dependence on mass 
also affects these fractions, with lower exponents predicting 
a higher fraction of stars with accelerations above a particu- 
lar value. Lower exponents correspond to a distribution with 
a larger number of massive objects. 

The top part of Table 2 shows distributions normal- 
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Figure 4. Probability contours sliowing tlie most likely mass 
and semi-major axis of a planet that could cause the acceler- 
ation 4.6±0.8 m s~^ yr~^ observed for GL 849 llButler et al.l 
|2006|) . The contours have been estimated using the planet and 
semi-major mass distribution and limits listed in the first line of 
Table [2] Contours are shown with probabilities, dP in log bins 
at 0.0003, 0.001 0.003, and 0.01. Flat eccentricity distributions 
were assumed, a) For an eccentricity distribution with maximum 
Bmax = 0.99. b) For emax = 0.3. The probability distribution 
is broader if the eccentricity distribution is broad. Objects with 
masses and semi-major axes with Mp/a? 0.03J\/,//(AU)^ arc 
most likely. 



ized and consistent with the known population of extra 
solar planets within about 5AU. None of the fractions of 
objects with accelerations above 10 and 25 are predicted 
high enough to match the fractions 23% and 12% esti- 
mated from the radial velocity surveys (|Cumming et al.l 
l2008l : lNidever et al.|[20o3 ). 

We consider the possibility that the known distribution 
of companions based on brown dwarf companion imaging 
surveys could account for the accelerations. In the set set 
of lines shown in Table [2] we show fractions of stars with 
accelerations above 5, 10 and 25 m s~^ yr~^ but normal- 
ized to_a^£ecentcompanion_^^ young clus- 
ter (|Kraus et al.ll2008l ). lKraus et all l|2008l ) find that 35% of 



stars have companions in the range 6-435 AU in the mass 
range ~ 10-500 Mj. We assume a stellar companion semi- 
major axis distribution consiste nt with a flat distr ibution 
in log a or an exponent /? = —1 (|Heacoxlll99 9: Kra us et al.l 
2008). For the companion mass distribution we use a shal- 
lower exp onent a in the range —0 . 5 to consistent with the 
study by iGrether fc Lineweaverl (|2006l ). The bottom part 
of Table [2] shows that fractions of objects with accelera- 
tions above 5, 10 and 25 m s~^ yr~^ predicted based on 
the binary companion distribution are also not high enough 
to match the ob served estimated fractions of 23% and 12% 
(|Cumming et a l. 2008; Nidever et al]|2002r ). 

We note that imaging surveys are necessarily incom- 
plete, particularly at small radii. Consequently a higher frac- 
tion of stars could host brown dwarf and stellar companions. 
However the fractions shown in the third set of simulations 
in Table [2] would have to be about 3 times higher to match 
the estimated fractions from radial velocity surveys. 

Both normalizing to the population of objects found 
from radial surveys and from imaging companion surveys 
predict a deficit of stellar accelerations. To explain the frac- 
tion of objects with detected line of sight accelerations an 
additional population of massive companions in the mass 
range a few to 100 M j is required. The extrapolation by 
ICumming et al.l l|2008l ) to larger semi-major axis using the 
mass and semi-major axis distribution suggested that 20% 
of stars would have gas giant planets within 20 AU. However 
here we find that this extrapolation is surprisingly conserva- 
tive and probably under-predicts the number of stars with 
observed radial velocity variations. 

We ask how many planets would be required to pre- 
dict the estimated fraction of stars with accelerations above 
10 and 25 m s"^ yr"^? The third sets shown in Table [2] 
show the fractions /lo and /25 predicted for stars assuming 
that 100% host a companion between 1 Mj and 80 Mj (the 
highest mass that would be considered a brown dwarf). If 
we restrict the companion between 4 and 20 AU and use 
siz e and period distrib ution consistent with that estimated 
bv lMarcv et al.l l|2008l ) then /lo and /25 are predicted to be 
similar to those estimated from radial velocity surveys. The 
last line in Table [2] shows the fractions assuming that 3 mas- 
sive companions are present in 70% of stellar systems. The 
planet inclinations in each system are assumed to be the 
same in the Monte Carlo simulation. This distribution also 
predicts fractions sufficiently high to be similar to those es- 
timated from radial velocity surveys. These two distribution 
wo uld predict more br own dwarfs than found by the survey 
bv lKraus et al.l l|2008l ). Imaging surveys may rule out mas- 
sive companions at large distances, implying that an even 
larger population of gas giant planets must exist in the out- 
skirts of extrasolar systems. 



4 SUMMARY AND DISCUSSION 

In this paper we have considered the statistics of accelera- 
tions that are detected in radial velocity surveys assuming 
that they are due to massive companions with periods that 
exceed the time span of the observational survey. We find 
that the acceleration scales with the planet mass divided by 
the square of its radius from the host star. The acceleration 
is independent of the stellar mass. As have previous studies 
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(e.g., iPatel et al.ll2007l ). we conclude that detectable accel- 
erations imply that there are gas giants or brown dwarfs in 
the outskirts (5 to 20 AU semi-major axis range) of their 
planetary systems. 

We have run Monte Carlo simulations to predict the 
statistics of the line of sight component of the radial stel- 
lar acceleration. We find that a distribution of companion 
masses and semi-major axes causes a probability function 
for the line of sight component of the stellar acceleration in 
the form of a constant minus a term proportional to the log 
of the acceleration. 

Given an object with a measured acceleration, Monte 
Carlo simulations exhibit a maximally likely product of the 
planet mass divided by its semi-major axis. The distribution 
is broader if the companion eccentricity distribution is wide. 

We find that extrapolation of the planet mass and pe- 
riod distribution found by radial velocity surveys within ~ 5 
AU and companion imaging surveys predicts too few sys- 
tems with detectable accelerations. This suggests that a sig- 
nificant fraction of stars, ~ 100%, hosts a gas giant or brown 
dwarf with semi-major axis between 4 and 20AU. Alterna- 
tively 70% of systems could host a few gas giant planets or 
brown dwarfs in a similar semi-major axis range. A large 
number of planets in the outer parts of extrasolar systems 
may not be unexpected if rapid clearing of extra solar sys- 
tems is required to expl ain debris disks with central clearings 
(|Faber fc Quillenll2007D . If imaging surveys rule out distant 
brown dwarf companions then we would infer that a large 
population of gas giants must be present in the outskirts of 
extrasolar systems. 

We have attempted to match an estimated fraction of 
stars ~ 23% and 12% with line of sight accelerations com- 
ponent above 5 or 10 and 25 m s~ ^ yr~^ estimated from 
the surveys bv lCumming et all l|2008l) ; iNidever et al] l|2002l ) . 
These fractions have not yet been precisely measured or re- 
ported. Stars showing a significant fraction of a single period 
or acceleration variations (as studied by iPatel et al. I l2007l 
with full orbital analysis) could be separated from those 
with constant accelerations. We have not considered the role 
of jitter and uneven time between observations and other 
noise sources in affecting the measurements of stellar accel- 
erations. Possible additional information could be learned 
from combining astrometric and radial velocity information 
for long period variations. The study here could be refined 
and broadened to more tightly constrain the population of 
objects in the outer parts of extrasolar planetary systems. 
As stellar and brown dwarf companions are ruled out from 
imaging surveys the fraction of objects at low mass can be 
better predicted. Better understanding of the statistics of 
accelerations will probably be needed to interpret forthcom- 
ing radial velocity and astrometric surveys of larger numbers 
of stars. 
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Table 1. Fraction of single planet mass and semi-major axis systems that would be seen 
with accelerations above 5, 10 and 25 m s~^ yr~^ 



Mp{Mj) 


a(AU) 


^max 


/5 


ho 


/25 


10 


10 


0.0 


0.36 


0.22 


0.00 


10 


10 


0.5 


0.34 


0.18 


0.02 


10 


10 


0.99 


0.29 


0.15 


0.03 


10 


5 


0.0 


0.47 


0.44 


0.34 


10 


5 


0.5 


0.46 


0.43 


0.32 


10 


5 


0.99 


0.45 


0.40 


0.26 


5 


5 


0.5 


0.42 


0.35 


0.14 


5 


5 


0.99 


0.40 


0.30 


0.12 


1 


5 


0.5 


0.15 


0.02 


0.00 


1 


5 


0.99 


0.13 


0.04 


0.01 



The fractions /s , fig and /25 are the fraction of stars predicted with line of sight accelerations 
above 5, 10 and 25 m s~^ yr~^ respectively. These fractions were measured from Monte Carlo 
simulations that considered a distribution of systems each with 1 planet with planet mass 
Mp (in Jupiter masses) (as listed in column 1), semi-major axis a in AU (as listed in column 
2) and with an flat eccentricity distribution in the range to e^^ax (as listed in column 3). 
Orbit inclinations and mean anomalies are randomly chosen. 



Table 2. Fraction of systems with a distribution of planet masses and semi-major axes 
predicted with accelerations above 5, 10 and 25 m s~^ yr~^ 



a (3 








^max 


/s 


ho 


/25 


Normalized to radial velocity searches 












-1.1 -0.6 


3 


100 


0.1 


100 


0.99 


0.05 


0.04 


0.02 


-1.1 -0.6 


3 


250 


0.1 


100 


0.99 


0.05 


0.04 


0.02 


-1.1 -0.6 


5 


100 


0.1 


100 


0.99 


0.04 


0.03 


0.01 


-0.6 -0.6 


3 


100 


0.1 


100 


0.99 


0.19 


0.14 


0.09 


-1.1 -0.4 


3 


100 


0.1 


100 


0.99 


0.09 


0.06 


0.04 


-1.1 -0.6 


3 


100 


0.1 


10 


0.99 


0.02 


0.01 


0.00 


-1.1 -0.6 


3 


100 


0.1 


100 


0.3 


0.06 


0.04 


0.02 


Normalized to an 


imaging search for companions 










-0.5 -1.0 


3 


250 


0.1 


500 


0.99 


0.08 


0.07 


0.06 


-0.5 -1.0 


3 


250 


0.1 


500 


0.5 


0.09 


0.08 


0.06 


0.0 -1.0 


3 


250 


0.1 


500 


0.99 


0.07 


0.06 


0.05 


0.0 -1.0 


3 


250 


0.1 


500 


0.5 


0.07 


0.06 


0.05 


Assuming 100% of stars have 


one massive planet 


or brown dwarf between 4 and 20 AU 




-1.1 -0.6 


4 


20 


1 


80 


0.3 


0.26 


0.19 


0.10 


Assuming 70% of stars have 3 


massive 


planets or 


brown dwarfs between 4 and 30 AU 




-1.1 -0.6 


4 


30 


1 


80 


0.3 


0.26 


0.21 


0.14 



The fractions /s, fig and /25 are the fraction of stars predicted with line of sight acceler- 
ations above 5, 10 and 25 m yr~^ respectively. The exponents a and correspond to 
a companion distribution dN oc M^a^ dMp da de. Minimum and maximum distribution 
semi-major axes a^in and amax are in AU. Minimum and maximum distribution planet 
or brown dwarf mass ratios are Mp^rnin and Mp^max in Jupiter masses. The eccentricity 
distribution is assumed to be flat with minimum and maximum Gmzn — and Gmax ■ 
The top set has distributions normalized so that the number of systems with masses 
between 0.3 and 10 Mj and periods between 2 and 2000 days is consistent with the 10.5% 
estimated by ICumming et al] ||2008| ) . The second set has distribut ions normalized so that 
35% percent have companions in the mass ran ge ~ 12 to 500 M j l lKraus et al.l l2008^. and 
with semi-major axis in the range 10-430 AU l lKraus et al.|[20o3). The exponent (3 for this 
set is chosen from the approximately flat distribution in log a " l|Heacoxl [l999l: [ Kraus et al.l 
20081). The mass expo n ent a is chosen to be in the range consistent with that found by 
Grether &: Lineweaved l|2006h . For the top 2 sets of distributions we assume that each 
system has at most 1 companion. The third set of simulations is computed assuming that 
all systems have exactly one massive planet or brown dw arf with semi - major axis between 
4 and 20 AU. Exponents a and (3 are those estimated bv iMarcv et al.l llooi). For the last 
set we assume that that 70% of systems have 3 giant planets or brown dwarfs in the region 
between 4 and 30 AU. 



